Time-Dependent Quantum Weak Values: 
Decay Law for Post-Selected States 
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Weak measurements offer new insights into the behavior of quantum systems. Combined with 
post-selection, quantum mechanics predicts a range of new experimentally testable phenomena. 
In this paper I consider weak measurements performed on time-dependent pre- and post-selected 
ensembles, with emphasis on the decay of excited states. The results show that the standard 
exponential decay law is a limiting case of a more general law that depends on both the time of 
post-selection and the choice of final state. The generalized law is illustrated for two interesting 
choices of post-selection. 

PACS numbers: 03.65.-w, 03.65.Ca, 03.65.Ta 



I. BACKGROUND 

Standard quantum mechanics famously describes two 
varieties of time evolution: the normal unitary evolution 
of the wave function, and the abrupt non-unitary col- 
lapse occasioned by measurement, when the system is 
projected onto an eigenstate of the relevant observable. 
The latter process involves an irreversible disturbance of 
the quantum system which alters its further evolution. 
Von Neumann's analysis of the measurement process, in 
which the quantum system of interest is strongly cou- 
pled to an external measuring device, explicates this non- 
unitary projection and unavoidable disturbance [? ]. A 
significant modification to this simple picture emerges, 
however, if one considers a quantum measurement car- 
ried out on an ensemble of systems. 

Consider a collection of N non-interacting systems pre- 
pared initially in the product state 



l*(^^) = l^)l|V'>2|^)3...|V')jV 



(1) 



and a set of identical observables {^j}, where Aj acts 
on the jth member of the ensemble, with corresponding 
Hermitian operators Aj . The ensemble average operator 
may be defined as 



N 



iW = (l/iV)5]i, 



(2) 



from which it is obvious that the expectation value of 
the average (A^^^) is the same as the average of the ex- 
pectation values (Aj). Let us now restrict to the special 
case that the N systems are identical, and prepared in 
the same initial state \ip), which may be achieved, for 
example, by making a strong, projective measurement of 
the same observable on every member of a large ensem- 
ble, and retaining in a sub-ensemble only those members 
that satisfy the required initial condition. It may then 
readily be proved [? ] that in the limit N ^ oo 



iW|vi/W)_(i)|M,W) 



(3) 



where (A) is the expectation value for a single member 
of the ensemble (the running index has been suppressed 
because all members arc identical). Furthermore, a mea- 
surement of A^^' will not disturb the system, or any in- 
dividual member thereof, in the large N limit. The fore- 
going properties are supported by detailed calculations of 
the actual measurement process, for example, based on 
von Neumann's model of a measuring device coupled to 
the entire ensemble [? ]. Thus the measurement yields 
non-trivial new information about the system, namely 
(A), without (in the large N limit) disturbing the mea- 
sured system at all. This result, which is contrary to 
the usual wisdom, may be understood intuitively by not- 
ing that (A) enters into Eq.Q as an eigenvalue of A^^' 
(though it is not in general an eigenvalue of A) . It is well 
known that repeated measurement of an eigenstate does 
not disturb the measured system. 

The key to understanding the foregoing result is readily 
apparent from inspecting Eq.([2]). If the measuring device 
couples to the whole ensemble with fixed strength, then 
the coupling to any individual member of the ensemble 
is reduced by the pre-factor of the right hand side, 1/A^. 
When A^ is very large, the coupling to individual sys- 
tems is very weak, and in the limit N —> oo the coupling 
approaches zero. Weak measurements imply small dis- 
turbance on the measured system, and in the limit being 
discussed, the disturbance will be arbitrarily small. It 
can be explicitly verified [? ] that the probability of a 
measurement disturbing any member of the ensemble ap- 
proaches zero like l/N. Nevertheless, information about 
the average is acquired. 

One may now ask what would be the outcome of a 
strong measurement performed on a given member i of 
the ensemble after the weak measurement has been per- 
formed to disclose some information about i. Although 
each member of the ensemble is prepared in an identi- 
cal state, subsequent strong measurements on i will not. 



of course, generally yield the same eigenvalue; rather, 
eigenvalues will be distributed with relative probabili- 
ties according to the Born rule. The initial ensemble 
may therefore be split into sub-ensembles according to 
the outcomes of the final strong measurement, and we 
are free to focus on a specific sub-ensemble. In other 
words, in addition to pre-selecting an ensemble of N 
identically-prepared systems, we may also post-select a 
sub-ensemble of < iV systems satisfying a final condi- 
tion too. In this paper, I shall consider the statistics of 
weak measurements carried out on specific pre- and post- 
selected sub-ensembles at intermediate times, i.e. after 
the (strong) pre-selection process and before the (strong) 
post-selection process. 

If the initial and final states of the whole ensemble are 
denoted |$i) and \'^f) respectively, and the correspond- 
ing single-member states by \ijji)j and \tpf)jj then the ex- 
pectation value of A for the jth member of the ensemble 
in the initial state may be decomposed as follows: 



(*d4l*.) 



k 



l(*.|V^M>jf(V^/,fcb^|*,)/(^/,fc|j*.) 



(4) 

where the complete set of single-member final states 
{\tpf.k)j} has been inserted. The first term in the sum- 
mand is the probability that the jth member of the en- 
semble will be found, on performing a strong measure- 
ment on j, to be in an eigenstate k of the (individual 
member) observable A. This probability will, in the limit 
N ^ oo,he the fraction of the ensemble that satisfies the 
generic pre-selection criterion plus the restricted post- 
selection of the particular eigenstate k. It thus defines a 
pre- and post-selected sub-ensemble. Equation^ there- 
fore shows that (Aj) of a single member j may be ex- 
pressed as a sum of all sub-ensembles (i.e. a sum over all 
post-selections k) of the total ensemble, with the relative 
fractions of the ensemble weighted by the corresponding 
quantities 



(5) 
Once again, because all members of the ensemble are 
identical, the index j is superfluous and will be dropped. 
The quantity 



{ipf.k\A\tlji)/{tljf^k\ipi 



(6) 



is known as the weak value of the operator, evaluated in 
this case for the specific post-selected eigenstate k. More 
generally, we may define the weak value of an operator 
A as 



w = {final\A\initial) / {final\initial) 



(7) 



for generic pre-selected initial and post-selected final 
states \initial) and \final) respectively. Weak values 



for pre- and post-selected ensembles are not expecta- 
tion values; rather, they are components of expectation 
values. They may lie outside the range of eigenvalues, 
and may even be complex, although the probability of 
post-selection yielding highly unusual weak values is nor- 
mally very small [? ]. Weak values have been the sub- 
ject of a considerable theoretical and experimental liter- 
ature. Good summaries of the interpretational aspects 
have been provided by Aharonov and Vaidman [? ] , and 
Aharonov and Rohrlich [? ], and will not be repeated 
here. The subject of weak values has stimulated research 
in topics as diverse as quantum tunneling [? ] and super- 
oscillations [? ] . It has also led to new experimental tests 
of quantum mechanics [? ? ] . In this paper I shall fo- 
cus on a specific and important example of weak mea- 
surements and post-selection, namely, the decay of an 
unstable system. 

Consider an ensemble of quantum systems prepared 
at time ti in identical quantum states |V'i)j and a sub- 
ensemble that is post-selected in state |V'/) at time if. 
The quantum weak value of an observable with Hermitian 
operator A is, at time t, ti < t <tf, given by 



w 



{^f\Wit-tf)Au{t-t,)\i^, 



{^Pf\Wit-tf)u{t-t,)\^p,) ' 

where U is the evolution operator for the system 



U{t) = e-'"' 



(8) 



(9) 



and H is the Hamiltonian. Most applications of Eq.® 
do not involve time-dependent evolution in the inter- 
val ti < t < tf. In this paper I shall consider explicit 
time dependence. (A rather different approach to time- 
dependent weak values involving weak energies of evolu- 
tion has been considered by Parks [? ]). 

A very basic property of quantum systems is the ex- 
ponential decay law, where the probability that a given 
system in an excited state at ti remains in an excited 
state after a time t is 



Pit) 



-27(t-ti) 



(10) 



where the decay constant, 27, is determined by the in- 
teraction strength between the system of interest and a 
set of "receptor" systems to which it is coupled. For 
example, the system might be an excited atom and the 
receptor system a bath of simple harmonic oscillators rep- 
resenting modes of the electromagnetic field. 

Equation PH)) is interpreted as follows: if a measure- 
ment is made on an individual system at time t, then 
the probability of finding it in the excited state is given 
by P{t). The measurement is understood to be a strong, 
projective, measurement, following which the state of the 
system is "reset" to \tpi) if it is indeed found that decay 
has not occurred. The question now arises of what will 



be found if a weak measurement is performed at time t. 
In the absence of post-selection, one would expect to re- 
cover the same exponential decay law. However, suppose 
that a post-selection were made at time i/; for example, 
consider the sub-ensemble of systems that have definitely 
decayed at time tf . What would be the decay law replac- 
ing Eg. ([TO)) , satisfying the dual constraints 



P{U) = 1 



(11) 



U{t) 



e-*"*/2 



from which one immediately verifies unitarity 
and the evolution property 



(19) 



(20) 



P{tf) - ? 



(12) 



To derive this weak measurement post-selection decay 
law, I shall consider an excited two-level atom coupled 
to a large bath of other two-level atoms, initially in their 
ground states, and compute the weak expectation value 
of the projection operator onto the excited state of at 
time i, subject to constraints (llip and p2)l . 

It is helpful to illustrate the calculation by first con- 
sidering a simpler time-dependent system. 



II. TIME-DEPENDENT WEAK VALUES 

Consider an electron of charge e at rest in a magnetic 
field B. The interaction Hamiltonian is 



where 



iJ= /i.B, 



/i = —ehS/m 



(13) 



(14) 



U{ti~t2)U{t2-t3) = U{h-t3). 



(21) 



Suppose at initial time ti the spin points in the +x 
direction, i.e. \ipi) is the normalized eigenstate of the 
spin operator ax with eigenvalue -|-1: 



IV'.) = 



72(1 



(22) 



a.lV',) =+|V'^)• (23) 

The projection operator onto the eigenstate ([221) is 



Prr 



1 A 1 



(24) 



^+ " V2 VI 1 
The expectation value of Px+ at time t > ti is then 

(V',|P,+ (t)|V'^) = {^r\UHt ~ U)Px+U{t - u)\-^,) 

1 (25) 

= -[1 -I- cosci;(i — ti)]. 



S = -{(Tx.ay^az) 



(15) 



and ai are the Pauli spin matrices. Suppose for simplicity 
that B lies in the z direction. Then Eq. (fT3| reduces to 



where 



H = fvjJCJz^ 



LO = eBz/m. 



(16) 



(17) 



The evolution operator Q is easily computed by ex- 
panding the exponential, using the relation 



trf = 1, 



(18) 



and summing, to find 



The projection operator onto the -1-1 eigenstate of ay 



p _ 1 A -i 



(26) 



from which it follows that 



{'lP,\Py+{t)\'lP,) = {'^^\U\t - U)Py+U{t - U)\ij,) 

1 (27) 

= -j[^ - s\nu;{t - ti)]. 

From Egs. f^Sl) and P?)) one may deduce the well- 
known result that the spin direction precesses with fre- 
quency oj in the x-y plane. 

The foregoing results may now be compared with weak 
measurements made at time t, with initial state |'0i) given 
by Eq.([221), and various choices of post-selection. For 
example, with 



l^/> 



V2 (-V ' V2 (-V ' V2 (l 



(28) 



corresponding to the eigenstates along the +?/, —x and 
the +x axes respectively, it follows using Eqs.® and (fTOl) 
that at time t, U < t < tf, the weak values w{Px+) of 
the projection operator Px+ are 



cos 



uj{t~ U) 



Uj(tf — t) . Loitf — t) 

COS —^ + sm —-^ 



uj{tf - U) . u){tf - U 
COS -— sm -— 



1 sm[uj{2t -U~ tf)/2] 

2 2sin[uj{tf-U)/2] 

cos[uj(t - U)/2] cos[uj{tf - t)/2] 
cosuj{tf — ti)/2 



(29) 



(30) 



(31) 



respectively. The corresponding weak values for P^^ may 
be derived in from the fact that, for the same choices of 

IV'j) and IV'/), 



be the eigenstate that the system would reach at time tf 
under unitary evolution from the initial state \'>pi). The 
result follows from the general property that 



{iPf\Wit~tf)Au{t-u)\,p,) 



(37) 



when 



\^f)^U{tf-U)\A) 



(38) 



which follows immediately by using Eqs. ipO)) and (PT|) . 

Expressions ([29 |) -([31 ]) display some curious properties 
typical of weak values. The weak expectation value for 
P can lie outside the normal range of probabilities [0,1]. 
It can, for example, take on both negative values, or val- 
ues > 1, for appropriate choices of tf, ti and t. Special 
interest attaches to (pij) . which is the weak expectation 
value of the projection operator onto the eigenstate of 
spin in the -l-a; direction, when this state is chosen for 
both pre-selection and post-selection. Naively, one might 
have supposed that w{Px+) would be fixed at 1 through- 
out, but in fact it is time-dependent. By way of illus- 
tration, if we choose condition (|34p . corresponding to a 
quarter of a precession period, then 



j{Px+)+w{Px-) = l. 



(32) 



It is also straightforward to calculate weak values for Py± , 
ax and ay. 

An important check on PT]) is to set 



'w{Px+) = -[1 + sin uj{t - ti) + COS ujit - ti)], (39) 

which takes the required value of 1 at t = t^ and t — tf, 
but in between is time-dependent and > 1. 



u>(tf — ti) = 2mr, 



1,2,3,.. 



(33) 



corresponding to an integral number of precession cy- 
cles, for then \ipf) = \ipi), and post-selection is trivial. 
Thus, applying condition ([33]) reduces ((3T|) to Eg. ([25]) as 
required. Similarly, if 



Uj{tf - t,) 

then \iljf) corresponds to the 
reduces to 



TT 

2 



(34) 
-X eigenvalue, and ([3D]) 



1 



[1 — COSbj{t — ti) 



(35) 



which is easily shown to correspond to the strong expec- 
tation value 



{'il},\Px-{t)\^^) = {i^^\UHt - U)P.-U{t ~ ti)\i,,). (36) 

Thus the weak and strong expectation values coincide 
for all t in the interval ti < t < tf when \ipf) is chosen to 



III. DECAY LAW 

I now come to the case of the decay of an excited state 
by considering an initially excited two-level atom coupled 
to a bath of 27V other two-level atoms initially in their 
ground states. The atoms are labeled by the index n. It 
is convenient to choose the ground state energies of all 
the atoms to coincide and be set to zero, and the excited 
states En to satisfy the relation 



En -Eo = nAE, -N <n<N, 



(40) 



i.e. the excited states are equispaced and distributed 
symmetrically about the excited state of the reference 
atom, labeled by n = 0. For simplicity I assume that 
the reference atom is coupled equally to each atom of the 
bath, and the interaction is described by the real constant 
Hamiltonian H. 

The Schrodinger equation for this system is equivalent 
to the coupled differential equations 



do = —i 2, HUnE 
n 



lAEt 



(41) 



-iHaQC 



inAEt 



(42) 



where a„ is the ampHtude of the excited state of atom n, 
and here and henceforth I set ft = 1 . Equations ([^T|) and 
(|32) niay be solved exactly (i.e. without the use of per- 
turbation theory) with the method of Laplace transforms 
[? 1 . One finds 



similarly, after some labor, but are not needed for the 
purposes of this calculation. 

Equation ([8]) may be rewritten, using relation ([20| plus 



uHt)^u{-t) 



(49) 



as follows 



ao(t) = C~ 



s + ^ coth I 



AE; 



(43) 



The inverse transform C^^ may be evaluated explicitly 
for n = in the limiting case N -^ oo, AE -^ 0, H ^ 
such that 



AE 



(44) 



where 7 is a finite constant and use has been made of the 
summation relation (jAll) in Appendix A. Then 



ao(t) 



- ^-7(t-*i) 



(45) 



Note that Eq.(l45|) is an exact solution (in the limit (|44|) ) 
for the amplitude of qq (t) . It follows that 



\ao(.tr=e 



2 = p-2-t{t-ti) 



(46) 



coinciding with the usual exponential decay law Eq.Q 
for the probability of finding atom in the excited state 
on performing a strong, projective measurement at time 
t. It has been derived here without the use of phase 
randomization or perturbation approximation. The fore- 
going analysis is thus ideally suited to generalization to 
weak values. 

To apply Eq. ^ we need to calculate the evolution op- 
erator U{t) for the system described by Eqs. ipTjl and (|^ . 
Restricting to the situation where only one atom at a 
time is excited, the evolution operator U{t) for the rele- 
vant subspace of the full Hilbert space of states will be 
a (2^^^ -I- 1) X (2iV -|- 1) matrix, the components of which 
may be calculated from the solutions of (PT|) and 
From Ea. (|l5)) we have 



Uoo{t) 



-ft 



(47) 



in the limit AE — > 0. Using this limiting solution, 
Eqs. (|iT|) and (|i^ may be solved to give 



Uno{t) = iH 



-ft+inAEt 



7- 



^AE 



(48) 



which is also understood as to be taken in the limit 
AE -^ 0. The other components of Unm{t) may be solved 



{^f\Uitf-t)Au{t~U)\i^i 



(50) 



If the operator A is chosen to be the projection oper- 
ator P| onto the excited state of atom 0, then Eq. ([5D|) 
describes the weak value for finding the atom undecayed 
at time t, given the condition that it is prepared in the 
excited state at time ti and post-selected to have decayed 
at time tf . A possible choice of final state is 



IV'/) = IV-fe 



(51) 



where the fcth atom only in the bath is excited, and the 
reference atom is in the ground state. In a physically 
realistic scenario, this corresponds to the atom in the 
ground state and a photon of energy Ek = kAE having 
been emitted. 

we note that P| is 



To evaluate ([50|) for this choice, 
represented by a (2A^ + 1) x (2A^ - 
components except 



(-Pt)oo = 1. 



1) matrix with all 



(52) 



The state (f/'fcl is represented by a row vector with all 
components except the fcth, which is 1 (for normaliza- 
tion). Similarly \Tpi) is a column vector with the i = 
component 1 and all others 0. Then equation (|50p re- 
duces to 



w 



UkQ{tf-t)Uoo{t-U) 



Uko{tf — ti) 
Using (gll) and (gH]) we find 



(53) 



w 



-7(t-ti) 



1 



-l(tf-t)+i{Ek-Eo){tf-t) 



1 _ e-'>'(*/-*^)+'(^---^o)(*/-*i) 



(54) 



The fact that w is complex is a familiar feature of weak 
values [? ] ; both the real and imaginary parts have phys- 
ical interpretations. In the case that Ek = Eq, i.e. the 
"photon" energy coincides with the reference atom's exci- 
tation energy, then (j54p reduces to the simple expression 



w = e-^(*"*') 



-7(*/-t) 



1 - e"'>'(*/~*«) 



(55) 



Equation 
cay law ps 
that 



551) generalizes the familiar exponential de- 
to the case of weak measurements. Note 



w = 1 
= 



t = u 

t - tf 



(56) 



as required, given that the reference atom is known to 
be in the excited state at t ~ ti and to be in the ground 
state aXt — tf. Furthermore, as i/ -^ oo, ([55|) reduces to 
the standard exponential decay law (j45|) , which is easy to 
understand because the system is known to have decayed 
in the limit t ^ oo, so post-selection of the ground state 
is redundant. 

Equation ([SO]) may be evaluated for other choices of 
operator A and final states |'0/)- One case of interest is 
A = P^ and the fcth component of (V'fc | given by 



(V'/l* 



iH 



■y + ikAE' 



Vfc. 



(57) 



Inspection of Eq. p5|) shows that this is the state to which 
the system approaches asymptotically as t — > oo, being a 
superposition of all excitations of bath atoms. It corre- 
sponds to the emission by atom of all possible photon 
energies. If this state is post-selected at the finite time 
tf, one must replace (|53p by 



where the summation is taken over all bath atoms. Equa- 
tions (pOj) and ([CT|) are consistent only if the summand 
contains both positive and negative contributions. The 
appearance of negative weak values for a projection oper- 
ator is another unusual but familiar aspect of weak mea- 
surement theory, as remarked in section 2. A full discus- 
sion of this topic will be presented in a future paper. 



IV. DISCUSSION 

The exponential decay law is one of the most funda- 
mental results of quantum mechanics, with wide appli- 
cability to atomic and nuclear physics. Concealed in the 
textbook discussion, however, is the fact that the expo- 
nential law ^ is actually a restricted case of a more 
general result for quantum de-excitation. The standard 
exponential result assumes that the system is inspected 
via a strong, projective measurement. By considering 
weak measurements the exponential law can be gener- 
alized, for example to (|55|) and ([55)) . depending on the 
choice of post-selection. 

Post-selection of sub-ensembles of atoms that are in 
their ground states at time t is conceptually straightfor- 
ward, and amenable to experimental test. 



APPENDIX A 



w = Uoo{t-ti) 



oo 

E 



-7(t-*0 



AB^O 



1 



Ukojtf-t) 
7 + ikAE 

p-27(t/-t) " 



/E 



Ukojtf — tt) 

7 + ikAE 



l_^-2j(tf^ti) 



(58) 



which also satisfies (I55|) and the condition that it reduces 
to (|45| when tf ^ oo. To arrive at Eg. ([58)1 I have made 
use of Eqs.lgl]), and (|XT|) and ((X2|) in Appendix A. 
Finally, I note that if one chooses |V'/) — \tpi), then 



w{P^) = 



Uooitf-t)Uooit^t,) 
U{tf - ti) 



= 1 



(59) 



in contrast to the time-dependent result ()3ip in the spin- 
precession case. In spite of this result, Eq. ([59)) does not 
imply that weak measurements of the bath atoms will 
show them remaining inert in this case. If we denote the 
corresponding projection operators for the bath atoms by 
P|„, then 



w{P^n) ^ 0. 
However, unitarity plus Eg. ([55)1 implies 



E"'(^Tn)=0 



(60) 



(61) 



The right-hand side of Eg. (|55)l involves the summa- 
tions 



lim AE V 



AE^o ^^ 7^ + k^AE^ 7 



(Al) 



Hm AE V 



AkAEt _ 

- - L—lt 



AE^o ^—i 72 + k^AE^ 7 



= -e-'^* it > 0) (A2) 



Equation (|Aip follows immediately from the relation 



aeY. 



k— — c 



-^ ,„ ^ ^„ = - coth (^) (A3) 

•j^ + k^jAE^ 7 VA£;/ ^ ^ 



(see, for example [? ]). Equation (|A2[) may be proved as 
follows: 



- ^^kAEt ^ 2 ^ cosh AM AE 

^ 72 -t- k^AE^ ^ AE^ k^+ 72/A^^ 7^ 
_ TT cosh[(7r - AEt)'y/AE] AE 



7 sinh(7r7/Ai?) 



7 



TT , , , TTJ , AE 



7 



[coth(— — - cosh7t — sinh7t] 



(A4) 



The summation on the first line of Eq. (|A4p may be for helpful comments, and Jeffrey ToUaksen for extensive 
found, for example, in [? ]. Finally, taking the limit assistance with the background theory. 



AE -^ 0, Eq.([A2|) follows. 
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